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ABSTRACT 


An  approximation  is  described  permitting  a  relatively  simple  analytic 
expression  for  various  statistical  averages  of  the  sound-pressure  field 
far  from  a  source  in  a  weakly  fluctuating  medium.  In  the  low-frequency 
regime  this  approximation  will  permit  valid  calculations  to  larger  dis¬ 
tances  than  the  conventional  eikonal,  or  geometrical  optics  on  straight 
Paths,  or  perturbation-theory  formulae,  provided  that  (kL)2(fec/c)4/3<  1, 
where  f'c/c.  and  L  are  the  fractional  sound-speed  fluctuations  and  correla¬ 
tion  length  respectively,  and  2n/k  is  the  sound  wavelength.  The  source- 
to-receiver  range  over  which  the  method  is  valid  is  R  =  L(l/kL)2(c/£c)2. 

The  method  is  used  to  estimate  mean  square  fluctuations  of  the  sound  pres¬ 
sure  and  coherence  lengths  in  an  isotropic  infinite  ocean,  and  it  is  pointed 
out  that  its  use  in  the  more  realistic  situation  of  cylindrical  symmetry 
is  not  appreciably  more  difficult.  More  careful  calculations  of  these 
quantities  are  presented,  ns  an  illustrate..  of  the  use  of  the  method, 
for  the  example  when  the  correlation  function  of  the  sound-speed  fluctua¬ 
tions  is  exponential. 
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I  STATEMENT  OF  THE  PROBLEM 


This  report  is  motivated  by  the  desire  to  understand  the  statistical 
fluctuations  observed  in  the  long-distance  propagation  of  low-frequency 
sound  in  the  deep  ocean.  There  are  temporal  and  spatial  fluctuations  of 
the  sound  speed  in  the  ocean,  caused  by,  among  other  things,  the  presence 
of  internal  waves,  and  these  lead  to  corresponding  fluctuations  in  the 
received  signal  from  a  sound  source.  Suppose  a  series  of  pulses  is  emitted 
from  a  sound  source  and  picked  up  by  a  distant  receiver.  During  the  very 
short  travel  time  of  a  single  pulse,  the  sound  speed  in  the  ocean  is  es¬ 
sentially  fixed,  but  between  successive  pulses  it  c;  n  change.  Thus  the 
variations  in  the  sequence  of  received  signals  probe  the  variations  in 
the  sound  speed.  Measurements  have  been  made  under  varying  conditions 
of  the  statistical  properties  of  such  sequences  of  received  signals, 
and  there  have  also  been  speculations  on  the  statistical  properties  of 
the  fluctuations  in  sound  speed.2  The  question  is  whether  or  not  these 
properties  can  be  connected  theoretically. 

From  the  above  physical  problem  we  abstract  the  mathematical  problem 
that  will  be  discussed  in  this  report,  as  follows:  what  are  the  statis¬ 
tics  of  sound  propagation  in  an  infinite  gaussian  random  medium?  Suppose 
we  have  a  CW  point  source  of  sound  of  wavelength  >.  located  at  the  origin 
in  an  infinite  medium  with  sound  speed  c  +  6c(x),  where  c  is  the  mean 
sound  speed  and  where  the  spatial  variations  6c(x)  are  small:  6c(x)  «  c. 
Suppose  6c(x)  is  given,  or  more  precisely,  suppose  the  correlation  function 


* 

References  are 


listed  at  the  end  of  the  report. 
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(1.1) 


and  where  k  =  2n/X  is  the  wav3number  of  the  emitted  sound.  For  a  point 

source  of  strength  S,  the  pressure  must  satisfy  the  boundary  condition 
2  2-*  3  — 

that  (V  +  k  )p(x)  -  S  5  (x)  as  x  ->  0;  we  shall  usually  choose  S  =  1  so 
that  we  are  dealing  with  a  source  of  unit  strength. 

The  parameters  that  will  be  relevant  to  our  discussion  are  >,  or  k; 

R  =  |x|,  the  distance  from  source  to  receiver;  and  L,  the  sound-speed 
correlation  length  such  that  C(x  -  x  )  vanishes  if  |x  _  x  |  >  L.  For 

the  purposes  of  explicit  illustration,  we  shall  sometimes  use  the  example 

-  I  x !  /L  4  2 

C(x)  =  C  e  ,  with  C  -  4k  (6c/c)  . 

Typical  numerical  values  of  these  parameters  are:  ).  ~  10  m,  L  ~  100  m, 

-4 

R  ~  a  few  tens  to  a  few  thousand  km,  and  6c/c  ~  10  .  These  values  are 

not  out  of  line  with  what  one  might  expect  for  the  oceanic  sound-propagation 
problem  that  motivates  us. 
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At  first  glance,  since  the  wavelength  is  small,  one  might  be  tempted 
to  think  that  geometrical  optics  provides  a  viable  approximation  with 
which  to  calculate  the  sound-pressure  field.  However,  in  order  to  obtain 
explicit  answers,  one  must  use  geometrical  optics  on  straight-line  propa¬ 
gation  paths,  and,  as  explained  in  Appendix  A,  the  condition  for  the  va¬ 
lidity  of  this  approximation  is  R  <  L(c/Sc)  ,  which,  with  our  typical 
parameter  values,  requires  R  ^  50  km.  This  is  too  short  a  propagation 
distance  to  satisfy  us,  so  we  are  compelled  to  find  a  better  approximation. 
How  this  is  done  is  the  subject  of  this  report. 


II  STATEMENT  OF  CONCLUSIONS 


To  spare  the  uninterested  reader  the  pain  of  paging  through  a  lot 
f  arithmetic,  we  shall  present  our  conclusions  here: 


(1)  An  excellent  approximation  to  the  average  received  pressure 
is 


where 


(p(x)> 


1  i[k-Z(k)/2k]R 

4  Tit  C 


(~  ■ 

(2r)3 


C  (  k  -  q) 


2  2 
-q  + 


(2.  i: 


(2.2 


and  where  C  is  the  Fourier  transform  of  the  sound-speed  corre^ 
lation  function: 

e(q)  =  /d3  x  e-1**’*  c(x)  .  (2.3 

3  2 

This  expression  is  valid  provided  (kL)  (6c/c)  «  1,  indepen¬ 

dently  of  the  distance  from  source  to  receiver. 

Since  Z(k)  is  complex,  (p(x)>  is  exponentially  damped  in 
distance.  We  have 


(p(x))|  =  e  ^ / 4nlt 


(2.4 


where,  for  the  case  of  spherical  symmetry,  where  C(q)  = 


-  2X 
C (q  ) 


“  '  1<in  f-i 


d(cos  9)C[2k  (1-  cosB)]  . 


(2.1 


_  -F 

In  particular,  for  the  exponential  example  where  C(x)  -  C  e 
we  find  (when  kL  »  1) 


" 


HI*  «IR  t 

,  _  f  ■:  —  -  r 


'■ 


(2) 


V 

a  =  CL/lk“  . 

1  2 

Thus,  if  C  =  Ik  (tc/c)  ,  we  get 

2  2 
a  =  k“L  (6c/c) 


(2.6) 


(2.7) 


so  that  the  damping  distance  is 

2 


(2.8) 


With  our  typical  values  for  the  parameters,  this  distance  is 
3,000  km. 


As  outlined  in  Appendix  A,  geometrical  optics  with  straight- 
line  paths  is  a  valid  approximation  for  distances  R  <  R0>  where 

H  =  min[L(c/6c)2  3,  L(kL) *  J].  Provided  that  R  R  ~  d  XL) 

Q  x  w 

(c/5c) ^  ^  >  1,  an  improvement  over  geometrical  optics  valid  out 
to  greater  distances  R  •  R^>  where 


It 


>  ■  (rMf 


(2.9) 


is  the  supereikonal  approximation,  in  which  the  pressure  is 
given  by 


-  C d3q  iq*  x  C  j  _i 

p(x)  =  if  - -  5  /  d' 

J  (2-)  Jq 


k  2  2  v  .  , 

-i?[q  -k  -l(S»x)-ieJ 


(2.10) 


where 


KS.o  f\ 

J  (2-)  Jo 


-isq'x  -ips(l-s)q 
ds  e  £ 


(2.11) 


and  where  V(q)  is  the  Fourier  transform  of  Y(x): 


~  _  /  3-*  -  iq  •  x  in  i 

V(q)  =  /  d  x  e  V ( x )  .  U*  1 

For  our  typical  parameters,  the  conditions  of  validity  peri  it 
the  use  of  this  approximation  out  to  a  distance  R  ~  3000  km. 


r> 


Note  that  this  range  of  validity  is  the  same  as  the  damping 

-1 

distance  a 

(3)  While  this  is  a  result  in  closed  form,  and  one  that  can  there¬ 
fore  be  used  for  numerical  purposes,  it  is  convenient  to  make 
a  further  approximation  in  order  to  obtain  a  more  explicit  ex¬ 
pression.  Under  the  condition  that 

kL  1  (2. 13) 


a  valid  approximation  to  Eqs.  (2.10'  and  (2.11)  is 


p(x)  =  (p(x) )  e 


X(x)-1/2(X  (x)) 


where 


X(x) 


-±f, 

4vJ 


3- 

d  y 


X  ik(  x-y  +  y  -  X  ).. 

— —  e  v 


x-y | |y  | 


(2.11) 


(y)  .  (2.15) 


(4)  From  Conclusions  (3)  above,  one  can  obtain  expressions  for  the 
statistical  averages  of  interest.  First, 


P(X1)P(X2)  [C(0)/16k]  [  iC(0)/8nk] (C+log  4kR) 

1  =  £3  £ 

<P(\))<P(^2))  (2.16) 

where  C  =  0.577...  is  Euler’s  constant.  Note  this  is  indepen¬ 
dent  of  x^-x^;  there  is  infinite  coherence  length  for  this 
quantity. 

Second,  for  the  case  of  spherical  symmetry,  the  mean  square 
pressure  is 


|P(x) l2* 

| <p(x))  |2 


exp  — 

4tt 


(2.17) 
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where 


C(q2)  =yV3xi1^C(x)  .  (2-18) 

To  illustrate  Eqs.  (2.16)  and  (2.17),  let  u',  use  our  exponen¬ 
tial  model 


C(x) 


C  e 


-  x 


l/L 


Then,  Eq.  (2.16)  becomes 


(p(x1)p(x,2)  >  i(CL3/k)(C+loK  4kR-in/2) 

<p(xi>N(p(x2)  ) 


(2.19) 


(2.20) 


while  Eq.  (2.17)  becomes 

(  [p(x)_[_J>  _  exp  RCL  .  (2.21) 

j (p(x) ) |2  k2 

If  we  set  C  =  4k4(£c/c)2,  these  expressions  are  replaced  by 


<P(X1)P(X2))  i.4(kL)3(f'C/c)2(C+log  4kR-in/2) 

-  =  e 

<fp(xi))(p(x2> 

and 

/|p(x)|2)  _ 4(kL) (5c/c)2(kR)  (2<2; 

|(p(x))|2 

One  must  keep  in  mind  the  conditions  of  validity  for  these  ex- 

2  2 

press  ions — namely,  R  «■  L(l/kL)  (c/6c)  •  Thus,  the  exponent 
in  Eq.  (2.23)  cannot  become  as  large  as  4.  Beyond  this  value 
of  R  the  exponential  growth  with  distance  ceases,  and  damping 
as  described  in  Conclusion  (1),  above,  sets  in.  As  a  result, 
Eq.  (2.23)  cannot  be  believed  at  distances  beyond  those  at 
which  the  exponent  is  of  order  one. 
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A 


if  dim, 


(5)  Finally,  while  the  remaining  quantities  we  set  out  to  compute 
are  more  difficult  to  obtain  explicitly,  rough  dimensional  es¬ 
timates  of  them  can  be  made. 


For  spherical  symmetry, 


P(X1)P  (X2)  :c<0)/(kL?][l-(x  /L)2Jlt 

X 

<p(x’1)  )(p(x^)  )* 


(2.21) 


where  x  =  x^  -  x^,  and  x^  is  the  component  of  x  perpendicular 
to  It  =  (x  +  x#^)/2.  For  cylindrical  symmetry,  if  It  lies  in  a 
"horizontal"  plane  transverse  to  the  symmetry  axis,  the  rough 
estimate,  Eq.  (2.24),  is  changed  to 


(p(xi)p*(x^) ; 

(P(x'1)  >(p(x2))1< 


Hr) 

k  L  L  '  V  '  '  H  ' 


(2.25) 


where  L  and  L  denote  the  horizontal  and  vertical  (trans- 
H  V 

verse  and  parallel  to  the  symmetry  axis,  which  we  call  "verti¬ 
cal")  sound  pressure  correlation  lengths,  and  x  and  x^  are  the 
corresponding  components  of  x. 

~  4  2  3  ~  4  2  2 

I f  we  estimate  C(0)  ~  k  (ftc/c)  L  and  C(0)  k  (fc/c)  L  L  for 

H  V 

the  spherical  and  cylindrical  cases  respectively,  we  find  for 
the  cutoff  distances  (i.e.,  the  pressure  coherence  distances): 


D  ~  )  J L/H  (c/6c) 


Dh  ~  >  7Lh/R  (c/6c)(Lh/Lv) 


DV  ~  >  7Lh/R  (c/6c)(Lv/Lh) 


With  our  standard  parameters,  all  of  these  are  of  order  10  m 
when  R  is  100  km. 
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Ill  THE  DERIVATION  OF  THE  APPROXIMATION 

A.  The  Perturbation  Series 

We  shall  now  show  how  to  derive  our  improved  expression  ior  sound 
propagation  in  a  random  medium.  The  fundamental  equation  for  the  propa¬ 
gation  of  sound  of  frequency  a  is 

(?“  +  k“)p(x)  =  V (x) p( x)  (3.1) 

where  k  =  jj/c,  c  is  the  mean  sound  velocity,  V(x)  =  2k  [6c(x)/cJ,  6c(x) 
is  the  (small)  difference  between  the  local  sound  velocity  and  c,  and 
p(x)  is  the  sound  pressure  field. 

If  V  =  0,  the  solution  of  Eq.  (3.1)  is 
ik|x| 

P  (x)  -  - -  (3.2) 

0  4P  |x| 

for  a  point  source  of  unit  strength  at  x  =  0.  Since  V  is  small,  it  makes 
sense  to  study  an  expansion  of  the  solution  in  powers  of  V.  To  that  end, 
we  introduce  the  Greens  function,  Gq(x  _  y) ,  satisfying 

(72  +  k2)GQ(x  -  y)  =  63(x  _  y) 

and  note  that  the  solution  of  Eq.  (3.1)  may  be  written  as 

p(x)  =  Pq(x)  +  J" Gq(x  _  x')  V (x) p(x  ,)  .  (3.3) 

Iterating  yields  the  desired  expansion  in  V: 

poo  =  P0U)  +  fi\  coa-xl)v<ii)po(;i) 

+  fax  dJx„  G  (x  -  x  )  V  (x  )  G  (x  -x  )V(x  )p  (x  ) 


A 


To  facilitate  further  manipulation  we  pass  to  momentum  space,  making 
use  of  the  Fourier  transforms 


r  -iqx  ,  v  1 

Go(ql=/,lxe  V*’  '"2 - 5 

^  k  -  q 

f  -iqx  1 

Po(q)  =  /  dx  e  P0(x)  =  — - - 

J  k  -  q 


(3.5) 


We  obtain 


p(q)  =  Y  P  (q)  ' 


P  (3)  =  f  II  tl3q  .V(q  ) 
n  ,1  1=1  i  1 


9  0  2  2  n  2  2 

(q“-k  )[(q  +  q  )  -k  ] . . . [ (q+ E q  .) -k  ] 


(n! ) J"  [J  d^q.Vfq  J 


n  n 

77  do .  Ml-Iii,) 

_ l _ i _ 

22  _22  n  2  2 

uQ(q  -k  )+o  t(q+q  )  -k  ]+  . .  ,+  q^  (q+  2^  )  -k  1 


The  last  step  makes  use  of  the  identity 


(3.6) 


fi(l-lQ() 


cl  •  •  •  cl 

1  n 


n  O'*  ^ 

L£o,i’ij 


(3.7) 


The  denominator  in  Eq.  (3.6)  may  be  rearranged  as 


D  =  (P  +Z  si  qi)2-k2  +  £  s^(l-s^)q3  +  X 


(3.8) 


where 


x  =  2  E  si(1"si)qi‘qi 


We  can  now  replace  the  integrations  over  the  o  in  Eq. 

tions  over  the  s..  We  have  only  to  remember  that  the 

0  and  1  and  are  ordered  s  <  s  .  Then 

i+ 1  i 


(3.6)  b\  integra- 

s  range  between 
i 


(3.10) 


and  therefore 


p  ( x ) 
n 


iq-x 


p  (q) 

n 


r  ,3  iq-xV 

"3—  —  -is  q.'xn 

I!  d  q  V(q  )  n  i  i  /7ds 

i  i  i 

J  q  c  J 

'  2  2  _  2  1  n-f  1 

• 

q  -k  +  £  s  ( 1 — s  )  q  +  X 

L  i  i  i  J 

(3.11) 


So  far  we  have  made  no  approximations  and  it  is  of  course  impossible 
to  sum  the  series  exactly.  We  shall  therefore  neglect  the  quantity  X  and 
show  that  with  this  simplification  the  series  becomes  summable.  Later  on 
we  shall  discuss  the  conditions  under  which  this  approximation  is  valid. 


Having  dropped  X,  we  see  that  the  integrand  of  Eq.  (3.11)  is  symmetric 


in  all  the  s  so  that  in  the  integration  over  s^  we  may  drop  the  restric¬ 
tion  s  >  s  and  absorb  the  factor  of  nJ  All  s  then  run  independently 
i  i+ 1  i 

from  0  to  1  and  we  have 


:u 

d  ci 


(2  )' 


i  •  •  -is  q  • x 
—  V  ( q  )  e  i  'i 
3  i 


2  2  11  2  ' 
q“-k“  +  £  s  ( 


2  1  n+  1 


n 

f  d  {«i 

.  -is  q...\ 

77 

i-1  J 

/  risi  - 3 

r  1  ( 21") 

1  (q  . )  •  i  1 

1 

r  •)  <>  n 

2  ~ 

q“-k“  +  £  i 

^  (l-s.)q 

L 

i  1  1 

1 1 


/3— t  _»  —  or. 


2  2  n  2 

-ipLq  -k  +  £s  (1-s  )q  -iej 
e  i  i  i 


H  d'*q  -is.q'.-x 

U  ds  - -  e  i  i 

i  =  l  1 


(2r) 


/d'q  -iq-x  .  f 

- 3  3  7 

(2r)  J  0 

IX 


2  2 

- ip  (q  — k  -  ie )  r—>  1  n 

dp  e  £  —  ( IP ) 

n. 


12..  - -  2 

d  q  -  isq  •  x-  ip  s  ( 1-s )q 

ds  - -  V(q)  6 

0  (2r) 


so  that  finally  we  may  write 


p(x) 


,  f  d'q  iq-x  f 

=  l 


iq • x  f  jo  -iP(q“-k2-I(P,x)-ie) 

e  I  dP  e 


with 


.  „  ,  dJq  fxn  -  i  sq  •  x-iP  s  ( 1-s  )q* 

I (6 , x)  -  /  - -  V  (q  )  1  ds  € 


*  (2r)  70 


(2.12, 


(2. 12) 


This  expression  may  appear  somewhat  unwieldy,  but  it  is  at  least  a 
closed  form,  and  we  shall  shortly  see  that  it  can  be  considerably  simpli¬ 
fied. 

At  this  point  we  note  that  if  we  had  dropped  all  the  terms  of  order 

2 

q  in  Eq.  (2.11)  we  would  have  obtained,  instead  of  Eq .  (2.12),  the  result 


I(P,x) 


(2.1  1) 


and  hence,  instead  of  Eq.  (2.12) 


This  expression  is  just  the  usual  eikonal  approximation  answer  lor  p(x), 
and  we  expect  that  our  new  answer,  Eq .  (3.12),  which  we  shall  call  the 
supe -•eikonal  approximation,  represents  an  improvement  over  this.  The 
nex.  question  to  discuss  is, when  is  this  answer  an  improvement  and  by 
how  much? 


1).  Error  Estimates 

Our  concern  here  is  not  to  show  precisely  what  are  the  limits  ol' 
validity  of  our  expressions,  but  simply  to  argue  heuristically  that  they 
should  hold  over  a  much  wider  range  than  the  usual  eikonal  result,  at 
least  under  a  wide  range  of  conditions  of  interest.  To  this  end  we  re¬ 
turn  to  the  exact  expression  for  p  (x),  Eq.  (3.11),  and  expand  the  denomi- 

n 

nator  in  powers  of  the  terms  quadratic  in  the  q^,  keeping  only  the  zeroth 
and  first-order  terms: 


ti 

C  d3^  iq.x  V'  f  d  ^i  „  “  i s  i  q  1  ’  x  . 
p ( x )  =  / - —  e  Lip  - -  v ( q  )  e  tfs 

J  (2DJ  nJ  1=1  ( 2r>) 3  1 


)  2  2  n+ 1 

(  (q  -k  ) 


(n+i) : 


Zs  ( 1— s  )  q  . 
i  i  l 

2  2  n+2 

(q  -k  ) 


(n+1) 

2  2  n+2 

(q  -k  ) 


(3.16) 

Here  we  must  remember  that  the  s^  range  from  0  to  1  and  are  ordered 

s  s  .  We  shall  call  the  three  contributions  to  p(x)  corresponding 
i  i  +  1 

to  the  three  terms  in  curly  brackets  Pg(x),  P^x),  ontl  P.,(x): 
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The  manipulations  that  allow  us  to  express  p^  and  p^  in  these  relatively 
simple  forms  are  straightforward  out  tedious  and  will  not  be  detailed. 

Now  pQ  is  just  the  normal  eikonnl  expression  for  the  pressure  field, 
p  is  a  correction  to  this  of  the  type  that  we  have  kept  in  the  super- 
eikonal  approximation,  and  p is  a  correction  we  have  thrown  away.  There¬ 
fore,  the  criterion  for  the  normal  eikonnl  to  be  valid  is  that  p  /pQ  1» 
and  the  criterion  for  the  supereikonal  to  be  valid  is  that  p^/p^  c  1. 

In  the  applications  we  shall  consider,  V  is  sufficiently  small  that 

we  may  set  yj k2+ J" *ds  V(sx)  equal  to  k.  Also,  we  are  interested  in  the 

behavior  of  the  various  p^  only  for  large  R  =  |x|.  Thus  the  indicated 
operations  5/fik  may  be  replaced  by  multiplication  by  R/2k  as  far  as  the 
leading  asymptotic  behavior  goes.  This  leaves  us  with 


1  R 2  /■*  (1-s) 

—  — ►  —  I  ds  - 

R-*0  2  /  s 

0  8k  -'o 


V2V(sx)  -  — —  §  ; 
8k2 


(3.20) 


_2 

p  3  k 

1  R-®- 


1  sl  (1'S 

/  dsi /  ds2  ~ r 

1  R  */0  -/0  1 


( 1-s  ) 

—  7V(s  x)-7V(s  x) 
1  2 


So 


ds  (1"S)  72V(sx) 


R 

=  —  a  . 

3k 


(3.21) 


To  discuss  the  behavior  of  these  ratios  as  a  function  of  R  we  need  to 

2  2  _ 

know  the  expectation  of  ft  and  $  .  Since  the  potential,  V(x),  is  assumed 
to  be  a  gaussian  random  variable  satisfying  (V(x)V(y))  =  C(x-y),  we  have 


<*2>  = 


/.l  /*1  (l-si 

7  dsi  /  ds2  “7" 

•'O  •'0  1 


(1-sJ  (1-s  )  22  4  — 

- -  s2  a  [V*C(X)]- 

S2  1  2  X=(s1-s  )x 


(3.22) 
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and 


<ii2> 


r1  f 1  f1  ,  r1  ,  ll-si 

=  I  ds  I  ds_  I  ds  I  ds_'  - - — 

•'O  •'0  •'O  •'0  “  fal 


(1-s.  )  (1-s^) 


/  / 
s  s  s  s 
12  12 


V2C(p)72C(0) 

p  =  (s  -s  )x 
12 

a  =  (s^-s^)x 


+  77  C ( p) 7  7  C (0) 

1  J  1  J  Ip  =  <s  -8^)5? 

a  =  (s2-s2)x 


+77  C ( p ) 7  7  C(0) . 

1  j  1  j  ip  =  (srS;>x  ) 


a  =  (s^-s 


I 


(2.22) 


Rather  than  attempting  to  evaluate  these  integrals  precisely,  we 

observe  that  C(x)  is  roughly  characterized  by  an  overall  scale, 

2  2 

[2k  (6c/c)]  ,  and  a  correlation  length  L  that  describes  the  distance 
over  which  it  falls  to  zero.  Then  each  factor  of  7^  in  the  integrals 
may  be  replaced  by  1/L  and  in  estimating  what  remains  we  may  regard  C(x) 
as  being  equal  to  [2k  (6c/c)]~  for  |x|  <  L  and  zero  otherwise.  Th^s 
gives  the  approximate  results 


( 


and 


4  ,  .  '2 


2>-7Kt)(S)  • 


(2.2  1) 


(2.25) 


the  factors  of  L/R  arising  from  the  fact  that  c((s^-s2)k)  will  be  large 
only  for  (s^-s^  S  L/R. 

We  now  can  put  the  error  estimates  in  a  manageable  form: 


P  2 

*1  __  R 


P  2 

*0  8k 


•v  /  2  v  R 2  /  2  6c  \  /T  1  1  6c  /  R  \3  2 

y<t>,_(2k  >(3.I6) 


1G 


and 


The  most  important  feature  of  these  ratios  is  their  different  If  dependence 
P/Pq  grows  much  more  rapidly  with  R  than  does  p^/p^.  Consequently  the 
approximation  of  neglecting  terms  like  p^  (the  supereikonal  approximation) 
should  be  good  out  to  much  larger  values  of  R  than  the  approximation  of 
neglecting  both  p^  and  Pr)  (the  normal  eikonal),  provided  that  the  coeffi¬ 
cients  in  Eqs.  (3.26)  and  (3.27)  allow  a  region  of  R  where  p_^/p^  1  while 

P^/P0  ^  1.  To  ill  istrate  the  sort  of  distances  involved,  we  choose 
2  “4 

L  =  10~m,  X  =  10  m,  6c/c  =  10  (these  are  roughly  the  parameter  sizes 
we  are  interested  in  for  deep  ocean  sound  propagation).  Thus 


and 


tor  R 


15  km 


for  U  £  3000  km  . 


The  different  R  dependences  of  the  two  error  estimates  have  a  dramatic 
effect  on  the  distance  for  which  the  two  approximations  are  valid.  It 
seems  clear  that  while  the  normal  eikonal  is  probably  inadequate  for  dis¬ 
cussing  most  long-distance  oceanic  sound  propagation  problems,  the  snper- 
eikonal  approximation  should  be  much  more  useful. 


C .  Further  Simplifications 

The  supereikonal  thus  constitutes  an  approximate  expression  in 

closed  font,  for  the  pressure  with  a  domain  of  validity  that  will  exceed 

the  domain  of  validity  of  geometrical  optics  on  straight-line  paths  pro- 
2  4/3 

vided  that  (kL)  (6c/c)  <  1.  While  these  expressions  are  relatively 
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straightforward,  further  simplifications  are  possible  when  the  propagation 
distance  is  large  and  the  wavelength  is  short,  and  we  wish  to  describe 
these  next. 


First  note  that  the  integral  over  d3q  in  Eq.  (3.12)  can  be  evaluated 
explicitly  to  yield 


p(x)  = 


(-in) 1/2 


8n 


7 

•'a 


3 


dg_  i(Pk+x  /4f3+pl(3,x)  +  ie 
3/2 


(3.28) 


When  k  and  x  are  very  large,  this  integral  may  conveniently  be  evaluated 
using  the  method  of  stationary  phase.  Let  us  define 


f(3)  =  Pk“  +  x2/43  +  pI(P,x)  -  3/2  log  p  . 


(3.29) 


Then  we  find 


p(x) 


2  e 


if<30> 


8T!  ^"(P0) 


(3.30) 


where  Pq  is  the  point  at  which  f'  vanishes: 


f'(P  )  =  0  . 

Next  we  note  that,  so  long  as  6c/c  «  1,  pQ  is  very  nearly  given  by 

PQ  =  x/2k  , 


(3.31) 


2  2 

the  point  at  which  p  +  x  /4p  has  its  maximum.  Therefore 


(3.32) 


our  result  is 


ikx 

e 

p(x)  =  - -  exp 

4  tlx 


(3.33) 


We  could,  of  course,  use  this  result,  which  is  now  only  in  terms  of  the 
one  complicated  integral  defining  the  function  I,  as  our  final  expression 
for  the  pressure.  But  we  can  do  even  better  than  this,  as  follows. 
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Let  us  recall  that  the  supereikonal  approximation  makes  its  first 
error  in  second-order  in  V,  and  to  first  order  our  expression  for  the 
pressure  must  be  exact.  Thus,  p^(x),  defined  by 


1  sn2  -L  f)3/2 


(3.34) 


must  be  io  correct  first-order  expression  for  the  pressure.  But  this 
integral  can  also  be  evaluated  by  the  method  of  stationary  phase,  since 
we  are  interested  in  large  values  of  k  and  x.  If  we  define 


g(ft)  =  gk2  +  x2/4g  -  3/2  log  (3  +  log  gl(£,x) 


(3.35) 


then  we  can  write,  approximately, 


—  ig(pQ> 
-  J  2  e 

v  '  8d 


(3.36) 


where  now  g  (PQ)  -  For  small  Cc/c,  we  again  have  =  x/2k;  hence 


we  find 


P.  Cx)  =  ~— 
1  4ux 


U  , ;)'  . 

J_  2k  \2k  / 


(3. 37) 


Combining  this  with  Eq.  (3.33)  yields' 


p(x)  =  PQ(x)  e 


p^^  (x)  /p0(x) 


where 


(3.38) 


p  (x)  =  - 

4-nx 


(3.39) 


is  the  zeroth  order  (in  V)  pressure  and 


-  r  3  elklx-yl  elky 

Jiu)VdyT^Tv<f,TT 


(3.40) 
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is  the  first  order  (in  V)  pressure.  Let  us  define 


X(x)  =  p  (x)/p  (x) 


1  f  3—  x  ik( | x-y |+y-x)  — 

=  —  /d  y  — -  e  1  y'  V(y) 

4tt  J  x-y  y 


(3.41) 


p ( x)  =  pQ(x)  eX*x) 


(3.42) 


It  will  be  convenient  in  what  follows  to  study  the  quantity 
p(x)/(p(x)).  Since  we  assume  the  fluctuations  in  sound  speed  to  be 
gaussian,  we  have 


,  ,  1  /2<X(x) 

( p(x) )  =  PQ ( x)  e 


(3.43) 


p(x)  _  ^X(x)-l/2 <X(x) )2 
(p(x)>  6 


(3.44) 


Over  the  propagation  distances  for  which  this  approximation  is  valid, 
(p(x)^  does  not  differ  appreciably  from  p^(x).  For  longer  distances, 
it  does,  and  a  quite  accurate  approximate  expression  for  it  is  derived 
in  Appendix  B.  We  find 


r 

i v  k2-Z(k) x 


<p(x) >  = 


(3.45) 


where 


/3-h 

/  n 


d  q  C ( k-q ) 
(2n)3  k2-q2+ie 


(3.46) 


and  C(q)  is  the  Fourier  transform  of  the  correlation  function: 


_  -iq*x 
x  e  C(x) 


(3.47) 
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F 


C(x-y)  =  (v(x)V (y )  ) 


(3.  18) 


From  Eq.  (3.44),  various  statistical  averages  of  the  received  pres¬ 
sure  are  readily  derived.  For  example, 


(pCx^pCx^)  (X(x1)X(x2)> 


(p(xi)  Xp(x2)  > 


\p(x^)p* (x^)  )  <X(xi)X*(x‘2)  ) 

(p(x^) )(p(x2) )* 


(3.49) 


(3. 50) 


Our  problem  is  thus  reduced  to  studying  averages  of  the  function  X(x), 
and  this  we  shall  do  in  the  following  section. 


IV  APPLICATIONS  OF  THE  SUPEREIKONAL  FORMULA 


Implementation  of  the  simplified  form  of  the  supereikonal  formula, 
derived  in  the  previous  section,  requires  the  evaluation  of  the  integral 
denoted  X(x),  or  what  is  perhaps  more  relevant,  the  evaluation  of  statis¬ 
tical  averages  of  the  integral  such  as  (X(x  )X(x,)>  and  (X(x  )x*(x,  )>. 
Before  attempting  this  in  general,  however,  it  is  of  some  interest  to 
see  how  the  conventional  geometrical  optics  limit  is  obtained  from  this 
formula,  so  let  us  indulge  in  this  digression  first. 

We  recall  that 


X(x) 


--  f 

4TI  J 


d3; 


l*~y 


ik(  x-y  + | y  -  x  ) 


V(y) 


(1.1) 


Let  us  choose  the  point  of  observation  x  to  lie  on  the  z  axis,  so  that 
x  =  (0, 0,  R) .  For  very  short  wavelengths— that  is,  in  the  limit  k  - 
the  only  important  part  of  the  region  of  integration  in  Eq.  (4.1)  is  neai 
the  z  axis.  Hence  we  may  expand  V(y)  around  the  z  axis: 


V(x,y,z) 


«  V(0, 0, z)  +  x 

Ox 


x=y=0 


1  2  5  V 

fiy2 


+  7y 


ov 

+  y  Oy 


b2V 


1  2  02V 

+  —  x  - 

2  2 

x=y=0  6x 


xy 


x=y=0 


fix2y 


+  •  •  • 


x=y=0 

(4.2) 


x=y=0 


We  further  note  that 


and 


z 


+ 


2  2 
x  +  y 

2z 


•  •  • 


R-z  + 


2  2 
x  +y 


2 (R-z) 


+ 


(4.2) 


(4.4) 
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finally,  the  z  integration  can  be  restricted  to  the  range  z  =  0  to  z  =  R. 
With  these  approximations,  the  integral  in  Eq.  (4.1)  can  be  readily  eval¬ 
uated.  We  note  that  only  terms  even  in  x  and  y  survive;  thus  the  terms 
in  i)V  x,  *)V/ty,  and  i)  V/£)x6y  in  Eq.  (4.2)  may  be  dropped,  and  we  obtain 

x<*)  =  f  d* / ^  oxi>f  )2) 

•'O  •'CU  •'CXJ 


|  V|  x=y=sQ 


2/2  _0_V 
+  2  \  ix2| 


2  i!2V| 


|x=y-0 


*y2  / 1 

x=y=0  '  ’ 


(4.5) 


Carrying  out  the  dx  and  dy  integrals  yields 

i  rK  1  rR  9 

X(x)  =  ~  2k  /  V|  dZ  +  - 7  / 

Jo  ,x=y=0  4k2  J0 


7  (z-z  /R)  dz 

x=y=0 


(•1.(5) 


which  is  recognizable  as  the  usual  formula  of  geometrical  optics  on 
straight-line  paths,  for  both  phase  and  amplitude  of  the  received  signal 


at  the  point  x. 

Let  us  now  return  to  the  more  general  situation,  and  attempt  to  com¬ 
pute  the  statistical  average  (X (x^ )X*(x^) > .  We  have 

<X(x1)X(x2)*>  =  (±f  C(yiry2) 

Kl  lx2l  ik(lx1-y1l+|y1  l-l^  |) 

K-yx  I  |yx  I  ’T*2-y2||y2l  e 

~ 1 k ( I x2“y2 I + I y2 I ~ I X2 I ) ■ 


(J.7) 

We  are  interested  in  the  situation  where  the  correlation  length  L  as  well 
as  the  wavelength  are  both  small  compared  to  1^  |  and  jxj;  thus,  ly^y  | 

lxif2^  and  I xi~X2 1  <<:  lxi  2 1  *  We  are  therefore  invited  to  change  to 
center  of  mass"  and  "relative'  coordinates  in  Eq.  (4.7)  and  to  assume 


that  the  relative  coordinates  are  small  compared  to  the  center-of-mass 
coordinates.  Thus  we  define 


x  +  x 

r  1  2  -  -  - 

H  = - ,  X  =  X  -X 
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y  +  v 

•/  *1  y2 


,  y  -  yl-y2 


(  I.K) 


and  we  write,  approximately,  when  |li|  »  lx  , 


I*!  I  3  l»  +  fl^,{  +  ^  «*x 


(4.9) 


with  similar  expressions  for  |xj,  |y  |,  |y^  | ,  |  x  -y  |  and  jx^-yj.  Then 


Eq.  (4.7)  becomes 


<X(x1)X*(x2)>  ^  |  —  J  / i'V  Jd3  it  'c(y) 


(*r 


iktlt  ( x-y )  +  It  •  y-h •  x  J 


(4.10) 


where  we  define  R  -  It  -  R  .  We  define  the  Foui  ier  transform  of  the  cor¬ 
relation  function  by 


C(q)  =  y*d3q"e  lq  x  C(x)  ; 

-(£)*  /dV(i^) 


(1. 11) 


i k ( It  -R)  *x  ~ 


—  /  * 

C(k(R  -  It  '))  .  (4.12) 


At  this  point  let  us  observe  as  an  aside  that  a  formula  analogous  t< 
Eq.  (4.12)  can  evidently  be  written  for  the  average  <X(x  )X(x9>): 


!x(x  )X(x  )>  =  (j-\  fd\'  (-?-_) 
12  4tt  )  J  IrV/ 


2ik(lt  ,+R,-H)  ~ 

e  c(0) 


(4.13) 
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This  formula,  we  note,  involves  only  the  correlation  function  at  zero 
argument,  and  does  not  depend  on  the  separation  x  =  x  -xy  between  the 
two  observation  points.  Evidently  the  situation  here  is  much  simpler 
than  that  which  obtains  in  the  calculation  of  \X(x^)X  (x^)  •  1  he  inte¬ 

gral  in  K(|.  (4.13)  can  in  fact  be  evaluated  analytically.  We  first  note 
tha  t 

d'V  =  — —  It  'dR  '  it  dll  (-1.1-1) 

R 

and  hence  we  find 

~  R/3  &  *1  ik  ( y — It/L? ) 

-  v  rc(o)  c  r  e  , .  ,r, 

<X(x  )X(xy)>  =  — —  J  dx  /  dy  - — -  (l.lf») 

-R/2  *8/2  V  "x 

where  we  have  defined  new  variables  of  integration  through 


R  '-U  " 

X  2 

Finally,  when  kR  »  1,  Eq. 

<X(x  )X(x  ) 
1  ^ 


u  '+r" 

and  y  =  — - —  • 

(4.15)  becomes 

)  =  lC(0-)-  (C  +  loK  4kR  -  ir/2) 
8nk 


(I.  1(5) 


(1.17) 


where  C  =  0.577...  is  Euler's  constant.  The  corresponding  pressure  aver¬ 
age  then  becomes 

p(xi)p(X2)  [  iC(0)/8rk]  (C+log  4kR-in/2) 

(p(x  )>(p(x ) > 

and,  we  recall, 

c(o)  =  yvx  C(x)  • 

There  is,  then,  a  distance-independent  amplitude  exp  C(0)/16k  that  is 
very  small,  and  of  no  interest,  and  a  small  phase  C(0)/8nk  (C+log  lkh) 
growing  logarithmically  with  distance,  of  almost  no  interest.  We 
recall  that  the  entire  approximation  is  valid  only  out  to  a  distance 


2  2  « 

K  ~  L(l/kL)  (c/6c)“,  so  the  growth  in  phase  in  limited  to  log  [ l/kL(e/tc)  1; 

furthermore,  as  we  remarked  earlier,  there  is  no  dependence  on  x  =  x  -x^. 

Let  us  now  return  to  Eq.  (4.12)  and  the  more  interesting — and  also 
more  difficult — average  (X(x^ )X* (x^)  ). 

To  begin  with,  we  shall  make  a  rough  estimate  of  the  integral  in 

Eq.  (4.12).  First,  suppose  the  correlation  function  C(x)  is  spherically 

symmetric;  then  C(q)  is  also  spherically  symmetric,  and  let  us  suppose  it 
2  2 

cuts  off  for  q  >  1/L  ,  where  L  is  the  correlation  length.  The  allowed 
region  of  integration  in  Eq.  (4.12)  is  then 

k2(2  -  2R  'it " )  £  1/L2 

A  f  *  // 

or,  if  0  is  the  angle  between  R  and  R  , 

0  £  1/kL  . 

If  kL  »  1,  this  is  a  small  angle.  The  integration  volume  is  therefore 
a  sausage-shaped  region  of  length  R  and  maximum  radius  of  order  R/kL, 
with  the  origir  and  the  point  R  at  opposite  ends  of  the  sausage. 

!  /  2 

Next  we  must  estimate  the  geometrical  factor  (R/R  R  )  .  This,  aver- 

2 

aged  over  the  region  of  integration,  turns  out  to  be  of  the  order  1/R  , 
as  one  might  expect  on  dimensional  grounds  alone.  On  the  basis  of  these 
comments,  let  us  first  estimate  (|x(R)|  );  that  is,  let  us  take  x  =  0. 

Then,  roughly,  Eq.  (4.12)  implies 


RC 

(kL)2 


4  3  2 

Since  C  ~  k  L  (6c/c)  ,  we  find 

< |X(R)  |2>  ~  k2RL  . 
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Next,  if  we  expand  the  exponential  in  Eq.  (4.12)  in  powers  of  x,  we  see 
that  the  linear  term  in  x  vanishes  by  symmetry,  and  the  quadratic  term  is 

2  j*  I  \  2 


(i;)  [kO?-R).xJ2f[k(H  -I? 


)J  . 


A  A  s\ 

Now  inside  the  sausage-shaped  column,  R  -R  is  nearly  perpendicular  to  R, 
and  is  very  small,  of  order  1/kL  i  1  length.  Thus,  only  the  component  of 

A 

x  perpendicular  to  R  survives,  and  we  may  estimate 


[ k(R'-R) • x] ~  ~  (x  /L)2  . 


As  a  result,  our  crude  estimate  gives 


(X(x^)X* (x^ ) )  ~  k2RL  -  Cx^/L)’ 


and  hence 


<P<X,)P*(I2»  ,  j  (W!  2  , 

~  "V  rl(t)  1i'<Vl>  1 1  • 

Thus  p(xl>p*(x2)s  vanishes  unless  k  RLCx^/L)  (6c/c)“  ^  1;  hence  the 
transverse  coherence  length  is 


D  ~  X 


-  (-) 

R  \  6c  / 


The  longitudinal  (i.e.,  parallel  to  R)  coherence  length  is,  in  contrast, 
infinite. 


We  may  also  estimate  (X(x^)X* (x^) )  in  cylindrical  rather  than  spher¬ 
ical  symmetry.  Let  us  suppose  that  C(q)  is  cylindrically  symmetric  in  a 

horizontal  plane,  and  cuts  off  when  q„  >  l/L  ,  and  q  >  1/L  ,  where  L 

H  H*  V  V  H 

and  Ly  are  horizontal  and  vertical  correlation  lengths,  respectively. 

Let  us  consider  the  case  that  it  lies  in  a  horizontal  plane.  Then  the 

region  of  integration  in  Eq.  (4.12)  is  a  flattened  sausage  of  length  R, 

3  2 

vertical  height  R/kL  ,  width  R/kL  ,  and  consequently  of  volume  R  /k  L  L  . 

v  ii  H  V 


The  factor  (R/R  R  )  still  averages  to  1/R*".  Then 


(X(xi)X*(x<>)>  ~ 


RC 


2 

k  L  L 
H  V 


where  xy  and  x^  are  vertical  and  horizontal  components  of  x  (the  longi¬ 
tudinal  component  of  x  does  not  enter,  as  before).  Next,  since  we  now 
have 


C 


‘1 

k 


2 

L  L 
11  V 


9 


we  find 


<p(xi )p*(x2> ) 
(pfx^ )  Hpfx^)  s* 


exp 


! *2«  L  (-)J 

/  xv  r 

/  XH  \2 

i  -  —  - 

(  —  1 

|  H\  c  / 

\  L  / 

L  v 

\  L 

H  i 

hence  the  vertical  and  horizontal  correlation  lengths  are 


Obviously  D^/D^  ~  E^/L^,  as  is  only  to  be  expected. 

These  crude  estimates  are,  perhaps,  useful,  but  it  would  obviously 
be  desirable  to  be  able  to  do  better.  Let  us  therefore  again  return  to 
Eq.  (4.12)  and  specialize  to  the  case  x  =  0.  Thus 


|x<5)|2>  =  (i)  /a'(fV)  £K-«'>) 


(1.  18) 


Let  us  suppose  we  have  spherical  symmetry.  Then  changes  of  variables 
just  like  those  used  in  obtaining  Eq.  (4.15)  yield  the  expression 


R/2 


< | X(R)  |  ) 


« fdx  f 

**  n  /o  •'rt  / 


dy 


-R/2  R/2 


2  2 

y  -x 


2k 


2  2 
2  y  -R  / 4 


2  2 

y  -x 


(4.19) 


where  we  have,  because  of  the  spherical  symmetry,  replaced  C(q)  by  C(q  ). 
The  distance  R  can  be  scaled  out  of  Eq.  (1.19)  to  obtain 


<|x(r)|2>  =  —  f y 

-i  *'i 


i 


2  2 

y  -x 


c 


2 

y  -i 

2  2 

y  -x 


) 


(4.20) 


n 

thus  < | X(R) j  )  is  simply  linear  in  R.*  A  further  change  of  variables 
finally  results  in 

<M>I2>  =  £ 


/  i  dZ  (sin"1  Vz)c(2k2z)  . 
/  v  z(l-z)  \  ' 

■/n 


(4.21) 


In  this  form  we  may  evidently  easily  compute  (|x|2>,  and  thereby  (|p|  » 

numerically  for  any  given  correlation  function  C.  Further  analytic 
progress,  however,  is  possible  only  for  sufficiently  simple  correlation 
functions. 


As  an  illustration,  suppose  we  look  at  the  example  of  an  exponential 
correlation  function: 


Then  we  find 


C(x)  =  C  e 


-x/L 


C(q)  =  8TTCL3 


2  2 
q  L 


(4.22) 


(4.23) 


This  statement  is,  of  course,  only  true  when  R  is  not  too  large.  The 
integrand  in  Eq.  (4.19)  should  really  be  modified  by  a  factor 

exp  -  4Ct(y  -  R/2) 

where  a  =  Im  Z(k)/2k  is  the  damping  produced  by  the  'self-energy  bubbles 
on  the  propagator  as  described  in  Appendix  B.  Thus  the  integrand  in  Eq. 
(4.20)  contains  a  factor 

exp  -  2QR(y-l)  ; 

hence,  if  uR  ^  1,  the  linear  dependence  on  R  becomes  modified.  However, 
when  QR  1,  R  is  so  large  that  the  supereikonal  approximation  fails,  and 
the  formalism  we  are  using  is  not  applicable  anyway.  Over  the  range  of 
validity  of  what  we  are  doing,  then,  i'  is  valid  to  neglect  this  factor, 
and  (|x(R)|2)  is  indeed  linear  in  R. 
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With  this  choice  of  C,  and  when  kL  is  large,  the  dominant  contribu¬ 


tion  to  the  integral  in  Eq.  (4.21)  comes  for  small  z;  hence  we  may  approx¬ 
imate  sin~^  s/~  in  Eq.  (4.21)  by  V* The  integral  may  then  be  evaluated 
analytically,  and  for  kL  large  we  find 


< |x(R) |2)  - 


RCL 

•> 

k“ 


4  2 

Setting  C  =  4k  (6c/c)  yields,  finally 


<  |X(H) | 2  >  =  4k2RL 


(4.21) 


(4.25) 


This  may  be  compared  with  the  very  crude  estimate  made  earlier;  we  see 
that  this  more  accurate  calculation,  for  the  exponential  correlation  func¬ 
tion,  differs  from  that  estimate  by  a  factor  of  4,  which  is,  we  note,  not 
totally  insignificant  numerically.  To  complete  this  example,  we  deduce 
that  the  mean  square  pressure  is 

( | p(x) | 2 )  _  [4( 6c/c)2k2L]R 

| (p(x)  }  j2 

2  s  2 

and  this,  we  recall,  is  valid  out  to  distances  R  ~  L(l/kL)  (i/6c)  . 

To  conclude,  let  us  summarize  what  we  have  found: 

(1)  For  an  arbitrary  correlation  function  C(x), 


<p(xi)p(x2)> 
p(xi) )(p(x2> ) 


exp 


C(0) 

16k 


exp 


1C(0) 

8TTk 


(C  +  log  4kR) 


where  C(0) 


/ 


d^x  C(x) , 


lx  +x„/2  . 


where  C  is  Euler's  constant,  and 


where  R 


(2)  For  an  exponential  correlation  function  C(x)  =  C  e 


-x/L 


<  ]p(iO  |2)  _  (CL/k^)R 
|\P(H)>|“ 


where  R  =  |x|  . 

(3)  For  an  arbitrary  spherically  symmetric  correlation  function 
C(x),  of  correlation  length  L, 


<p(x1)p*(x2)>  [5(0)/,kL,2Hl-(x/L)2»] 


<p(x^) )<p(x2) ) 


T.  \  ~  ® 


where  x  =  x  -  x# ,  and  x^  is  the  component  perpendicular  to 
x  +  S  K.  Thus,  the  transverse  coherence  length  is  of 
the  order 


D  ~  L 


kL 


ys 


C(0) 
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Appendix  A 


GEOMETRIC  OPTICS  WITH  STRAIGHT-LINE  PROPAGATION 


We  wish  to  briefly  review  the  conditions  under  which  geometrical 
optics  with  straight-line  paths  are  valid.  This  is  also  the  same  as 
the  result  obtained  from  the  supereikonal  formula  in  the  very-short- 
wavelength  limit,  as  outlined  in  Eqs.  (4.2)  through  (4.6).  The  first 
term  ii.  Eq.  (4.6)  is  what  is  obtained  if  all  quadratic  terms  in  the 
momenta  are  neglected  in  the  derivation  of  the  supereikonal  formula; 
this  is  what  is  known  as  the  conventional  eikonal  approximation. 


The  wave  equation, 


(V2  +  k2)p(x)  =  V(x)p(x) 

may  be  rewritten  as  the  two  equations 

-  -  _  2  2x  - 

-2k  x-76x  -  (76X)  -  —  *7  log  A/A 

x  0 

2  _  2 

+  7  log  A/A  +  (7  log  A/A  )  =  V 

0  0 

and 

2  x  A  — • 

7  6x  -  2  —  •  76x  +  2k  x-7  log  A/A 
x  0 

+  2V6x*7  log  A/Aq  =  0 


through  the  definition 

-i  log  p(x)  -  kx+6x(*)  +  i  log  4nx  -i  log  A/Aq 


(A-l) 


(A-2) 


(A-3) 


(A-4) 


where  Aq  =  l/4nx,  and  where  6x  and  A  are  real.  Equation  (A-4)  may  also 
be  written 
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ikx 

p(x)  =  —  •  (A/A  )  e  1  X 

4Tix  0 


(A-  5) 


Thus,  6\  is  the  change  in  phase  and  A/A^  is  the  change  in  amplitude,  d 
to  the  presence  of  the  potential  V(x). 


ue 


Now  in  Eq.  (A-3)  let  us  assume 


(a)  (76 x)  «  V 


(b)  7  log  A/Aq  «  V 


(c)  —  •  7  log  A/A  «  V 

x 


(d)  (7  log  a/Aq)  V 


and  in  Eq.  (A-4)  let  us  as 


sume 


—  —  2  x  — 

(e)  76x  *  7  log  A/A  «  V  6\»  —  •  V6x, 

0  x 


kx-7  log  A/Aq  . 


Then  Eqs.  (A-3)  and  (A-4)  reduce  to 


and 


-2k  x- 76x  =  V 


2  ^ 

7  6x  -  2  -  •  76 X  +  2k  x-7  log  A/A  =  0 
x  0 


(A- 6) 


(A-7) 


which  are  the  equations  of  geometrical  optics  with  straight-line  paths, 
or  of  the  conventional  eikonal  approximation.  Their  solution  is  trivially 
seen  to  be 


and 


jx(r,e,e*)  =  -  ±fdr'  V(r',6, 


«0 


(A-8 ) 


,  A(r,9,(rf)  1  r  /  /  /  r  '2  \  _2 

l0g  A~('r,e,«0  =  ~2f dr  (r  )7  V(R  ,e,«  .  (A-9) 


Given  these  solutions,  we  can  now  go  back  and  chec!  under  what  conditions 
Assumptions  a  through  e  were  valid. 
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First  we  note,  from  Eq.  (A-8) ,  that 


.  v 

or  k 


but  that 


_  .  R  V 

7  6X  ~  7  r 

x  L  k 


where  7  means  the  gradient  transverse  to  the  path  of  integration.  Hence 


„2,  R  V 

v  6x~- -  . 

L 


Similarly,  we  see  from  (A-9)  that 


6  ,  R  V 

—  log  A/A 
Or  0 


but  that 


and  that 


2 

2 

L 

k 

2 

R 

V 

3 

2 

L 

k 

2 

R 

V 

4 

2 

L 

k 

Using  these  estimates,  our  assumptions  a  through  e  now  require  that 

v2 


/  R  Vf 

“  (li) 


«  V 


2 

R  V 


<b>  -7  T  v 


4  2 

L  k 


(c)  -  •  —  —  «  V 

R  L2  k2 


(d) 


,  2  3 

ftj) 


«  V 


L  k 


(e) 


R  V  R V  V  R 

L  k  '  ,3  ,2  k  2 

L  k  L 
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Among  these  conditions  we  note  that  a  and  b  imply  d;  that  a  implies 
e  and  that  c  simply  requires  kL  <"*  1.  Thus  only  the  two  conditions  a  and 
b  count,  and  they  require 

(a  ')  R  «  (kL)  v  T Tv 

(b 7)  R  «  L(kL)  . 

These  yield  very  limited  ranges  of  validity  indeed;  however,  we  may  im¬ 
prove  them  by  recalling  that  we  are  really  only  interested  in  statistical 
fluctuations  in  the  pressure.  Thus  we  are  justified  in  multiplying  the 
left-hand  sides  of  Conditions  a  and  b  by  1/  where  N  ~  R/L  is  the  num¬ 

ber  of  traversed  inhomogeneities.  Then  a  and  b  become 


(a  ) 


L  k 


which  yield  the  limits 


and 

respectively. 

To  conclude, 


(a")  R  «  L(k2/V)2/3 
(b")  R  <<  L(kL)43 

2 

let  us  recall  that  V  =  2k  6c/c;  thus  we  have  to  have 


(a'  )  R  «  L(c/6c)2//3 

■1/T 

(b"j  R  «  LC2TTL/}  ) 


in  order  to  be  able  to  use  geometrical  optics,  or  the  eikonal  approxima¬ 
tion,  or  whatever  one  wants  to  call  it. 

Condition  a",  we  note,  was  also  obtained  in  our  derivation  of  the 
supereikonal  approximation,  os  the  condition  for  the  validity  of  the 
ordinary  eikonal  approximation. 


36 


Finally,  for  our  standard  values,  the  validity  conditions  of  Reo- 


metrical  optics  state  that  we  must  have 


(a") 

2 

R  «  10  • 

(ioV/3 

(b") 

2 

R  «  10  • 

<60)4/3m 

m 


-  50  km 

25  km. 
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Appendix  B 

CALCULATION  OF  SOUND  PRESSURE  AT  GREAT  DISTANCES  FROM  A  SOURCE 
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Appendix  B 


CALCULATION  OF  SOUND  PRESSURE  AT  GREAT  DISTANCES  FROM  A  SOIRCE 

\Ve  wish  here  to  briefly  outline  a  technique  different  from  that  dis¬ 
cussed  in  the  text  for  handling  calculations  of  the  sound  pressure  at  very 
great  distances  from  a  source.  This  technique  permits  a  very  good  calcu¬ 
lation  of  <p(x)>,  but  is  not  well  suited  to  deriving  useful  analytic  re¬ 
sults  for  more  complicated  averages  such  as  (|p(x)|  • 

The  idea  is  based  on  the  use  of  the  methods  of  quantum  field  theory, 
and  it  has  appeared,  albeit  in  what  seems  to  be  unnecessarily  complicated 

language,  in  the  literature.4 

We  start  with  the  perturbation  series  for  the  pressure,  Eq.  (3.4) 
of  the  text: 

p(x>  =  P0<x>  *  f*\  00(x:.x1)V(xi)p0(x1) 

♦  Ji\  /<<%vxVv<VVvVnVVV 


+  . 


(B-l) 


and  we  note  that,  for  unit  source  strength, 


PqU)  =  Gq(x)  . 


(B-2) 


Let  us  calculate  the  average  pressure  from  Eq.  (B-l).  Recall  that  \  (x) 
=  0,  and  let  us  assume  that  only  two  point  correlations  matter,  so  that 


<V(x  ) . . . V(x  ))  =  0  if  n  is  odd 
1  a 
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while 


\V( x  ) . . . V(x  ) ) 
1  n 


2^J  C(x-x)..  .  C  (  x  -x  ) 

1  2  li-l  n 


pc miu tut  ions 


if  n  is  even. 


Tlius  wo  find 


p ( x >  -  Go(x)  +  fd\  f'%  Go(x-xi)G0(xrx1,) 


C(X-X2)C0(X2) 


(15-3) 


This  is  obviously  precisely  the  perturbation  expansion  of  the  propagator 
for  a  particle  whose  free  propagator  is  given  by  G^fx),  which  can  emit 
and  absorb  another  particle  whose  propagator  is  C(x).  The  set  of  Feynman 
diugiams  comprising  p(x)  is  shown  in  Figure  II—  1,  and  the  Feynman  rules 
for  calculating  the  contribution  of  any  diagram  are  evident  from  Kq.  (15-2). 


<p(xl>  = 


+  +  r&Zb  +  ... 

FIGURE  B-1  FEYNMAN  DIAGRAMS  FOR  THE  "PROPAGATOR"  <p(x)>.  The  straight  line 
represents  the  propagator  Gq;  the  wavy  line  represents  the  propagator  C. 


Similar  expressions  can  be  easily  derived  f i  om  Kq.  (15-1)  for  other 
averages  of  interest.  For  example,  the  quantity  (p(x^)p(x  ) )  is  repre¬ 
sented  by  the  diagrams  shown  in  Figure  15-2,  and  is  evidently  a  vertex 
function,  in  which  two  "particles”  emitted  at  the  origin  propagate  to 
the  points  x^  and  x^  emitting  and  reabsorbing  the  "c  type"  particles  as 
they  go.  The  quantity  \p(x^ )  p(x^)’*  >  is  represented  by  the  same  set  of 
Feynman  graphs,  except  that  one  or  the  solid  lines  now  stands  for  propa¬ 
gators  G*  rather  than  Gq. 
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0 


1 


<p  (xi  )  p  (x^)  >  — 


FIGURE  8-2  FEYNMAN  DIAGRAMS  FOR  THE  "VERTEX  FUNCTION"  <p ( x  1 ) p < x 2 P 


Now  let  us  return  to  Eq .  (H-3).  Since  6c/c  is  so  small,  the  quan¬ 
tity  C  is  very  small.  Hence  the  dominant  set  of  graphs  contributing  to 

■  4 

p(x)  will  be  those  shown  in  Figure  H-3.  These  are  trivially  summed  if 


we  go  to  momentum  space: 


<P(x)>  = 

J  (2n)J 


X  _ 1_ 

2  2 
q  -k  +J 


1  P  t<3q  '  C  (q*-q* ')  + 
q2-k2+ie  (2n)3  q'2-k2+ie 


d  q  iq  •  x  1 

- —  p  M  - - 

3  2  2 

(2”)  q~-k  +  £(q)  +  ie 


( B— I  > 


where 


(  2  tt)  k  -q 


(H-5) 


<p(xl>  » 


FIGURE  8-3  DOMINANT  FEYNMAN  DIAGRAMS  CONTRIBUTING  TO  <p ( x )>  WHEN  ftc/c  IS  SMALL 


11 


At  very  large  distances  It  -  |x|, 
from  the  singularity  in  q  in  the 


the  leading  term  in  (p(x)  will  come 
Fourier  transform;  this  occurs  when 


q2-k“  +  I(q)  =  0 

which,  When  I  is  small,  is  approximately  given  by 

•>  2 

q“  =  k  -  I(k)  . 


(B-6) 


(B-7) 


Thus  we  arrive  at  the  answer: 


(p(x)>  = 


i  •/k“-I(k)U 


4rll 


plkH  -Ul(k)/2k]R 

-  -  p 

4  nit 


(B-8) 


In  view  of  the  ie  in  the  denominator  of  the  integrand  in  Eq.  (B— 5) ,  Z(k) 
has  an  imaginary  part  as  well  as  a  real  part.  This  gives  rise  to  a  damp¬ 
ing  of  ;P(x)>.  If  we  forget  the  uninteresting  small  phase  charge  produced 

by  He  £(k),  we  can  write 

iklt 


_  ,  e  -qk 

<P(x)>  =  e 


(B-9) 


where 


o  =  -Im  Z(  k)  =  Y1 


f- 


—  C(k-q ')  Mq  '2-k2) 


(2r) 


(B-10) 


For  example,  if  we  choose 


C (x)  =  C  g 


-  x|/L 


(B-ll) 


we  find,  from  Eq.  (B-10), 


a  =  CL/4k 


(B-12) 


The  use  of  the  Feynman  diagram  technique  is  easy  and  valid  for  cal¬ 
culating  (p(x)>.  It  turns  out  that  it  is  also  easy  to  calculate 
(p(xi>p(x2> 


(p(xi> >{p(x2)> 


;  the  lowest-order  diagram  shown  in  Figure  B-4  is  all  that 
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<p(x1)  p(x2)  > 
<p(x1)><p(x2)> 


<  *  <: 


FIGURE  B-4  THE  ONLY  DIAGRAMS  NEEDED  TO  CALCULATE  <p(x 
WHEN  iSc/c  IS  SMALL 


<p(x1  )p(x2))/<p(71  |><pU2)> 


it  is  necessary  to  include.  The  result  obtained  by  doing  this  is 


<p(xi)p(x2) ) 

7('xi)\P(x2)  =  1  +  <X(^)X^2)> 


(B-13) 


where  X(x)  is  as  defined  in  Eq.  (3.41).  Thus,  the  Feynman  diagram  tech¬ 
nique,  for  this  quantity,  coincides  exactly  with  what  we  obtain  from  the 
supereikonal.  (We  recall  that  the  expansion  exp  <X( x  )X(x  ) ^ 1 +  (X(x  )X(x 
is  quite  valid  for  this  case.) 

It  is,  however,  not  easy  to  calculate  (p^p^)*)  by  this  method. 
The  first  Feynman  diagram  coincides  with  the  first  term  in  (X(x  )X*(x  )' 
of  the  supereikonal,  but  here  it  is  not  valid  to  stop  after  one^rm/even 
when  bc/c  is  small.  An  infinite  series  of  diagrams  is  needed,  and  (a  part) 
of  such  a  series  is  what  is  provided  by  the  supereikonal  formula. 
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